Nuclear isotope thermometry 
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We discuss different aspects which could influence temperatures deduced from experimental isotopic 
yields in the multifragmentation process. It is shown that fluctuations due to the finite size of the 
system and distortions due to the decay of hot primary fragments conspire to blur the temperature 
determination in multifragmentation reactions. These facts suggest that caloric curves obtained 
through isotope thermometers, which were taken as evidence for a first-order phase transition in 
nuclear matter, should be investigated very carefully. 



I. INTRODUCTION 



Due to the short range attraction between nucleons, 
nuclear matter is a Fermi liquid jl[ at low temperature 
and is expected to undergo a phase transition to a nu- 
cleonic gas within a mixed phase region bounded by a 
critical temperature of order 15 MeV [[||||. Experimen- 
tal investigations of this phase transition have focused on 
a variety of experimental observables ranging from the 
mass, charge or multiplicity distributions for the emitted 
fragments QJjj] to observables sensitive to the tempera- 
ture of the system g^]. 

Temperature measurements, in particular, have been 
performed to search for evidence of the enhanced heat 
capacity predicted by statistical model calculations re- 
flecting the latent heat for transforming the Fermi liquid 
to the nucleonic vapor |^-^]. For example, the Statistical 
Multifragmentation Model (SMM) || predicts a plateau 
of roughly constant temperature of T w 5MeV for exci- 
tation energies of E* /Am 3 — 7 MeV. At these excitation 
energies the model predicts a mixed phase consisting of 
fragments (liquid) and nucleons and light particles (gas) 
corresponding to a mixed phase equilibrium. This is fol- 
lowed at higher excitation energies by a linear rise in the 
temperature with excitation energy as expected for a gas 
of small nuclei having negligible internal heat capacity 
||. Similar effects are predicted by the Microcanonical 
Metropolis Monte Carlo (MMMC) model JhJ. 

This trend was qualitatively re pro duced in some ex- 
periments 0, but not in others |l l||l4| . An essential 
part of these measurements is the determination of the 
temperature of the fragmenting system. Temperatures 
were extracted from the isotopic abundances of helium 
and lithium fragments, using the isotope thermometry 
method proposed by Albergo et al. yl| . The idea of the 
method is to determine the double ratios of the yields 
of four suitably chosen isotopes, (Ai,Zi), {A\ + l,Zi), 
(A 2 ,Z 2 ),(A 2 + 1,Z 2 ), 



where the Y are the yields of the different isotopes, C 
is a constant related to spin values and kinematic fac- 
tors, AB = B{Ai,Zx) ~ B(Ai + \,Z X ) - B(A 2 ,Z 2 ) + 
B(A 2 + l,Z 2 ) is obtained from the binding energies of 
the isotopes appearing in Eq. ([!]), and Ti SO stands for 
the temperature deduced from this isotopic thermome- 
ter. In the case of the He-Li thermometer employed in 
ref. §, A x = 6,Zi = 3, A 2 = 3,Z 2 = 2. For the C-Li 
thermometer, more recently considered by Xi et al. [ [l4|| , 
A x = 6, Zi = 3 7 A 2 = 11, Z 2 = 6. For the Carbon ther- 
mometer studied in this work, yli = 12, Z\ = 6, A 2 — 
U,Z 2 = 6. 

However, there are a few aspects which should be care- 
fully analyzed when one wants to compare information on 
the breakup configuration of an excited system formed in 
a heavy-ion collision to multifragmentation models like 
the SMM approach. Some of these points are addressed 
below. In sect, [n] we briefly discuss the assumptions un- 
derlying this method. Variations in the temperature of 
the breakup stage, where the hot primary fragments de- 
couple from the system, are intrinsic to finite systems 
and are explored within the SMM approach in sect. III. 



An analytic description of temperature variations is de- 
veloped in the grand canonical limit in sect. IV; this de- 



scription is consistent with the results from the SMM. 
In addition, there are finite size effects, discussed in sect. 
0, that make the concept of an overall chemical potential 
somewhat inaccurate. The influence of secondary decay 
is discussed in sect. VL Conclusions are drawn in sect. 



VII 



y(Ai,Zi)/y(Ai + i,Zi) 

Y(A 2 ,Z 2 )/Y(A 2 + 1.Z 2 ) 



Cexp(A J B/T lso ) (1) 



II. UNDERLYING ASSUMPTIONS 



The basic physical hypotheses of the isotope thermom- 
etry method are: 

1. an equilibrated source is formed after the most vi- 
olent stages of the reaction and it decays simulta- 
neously and statistically, 
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2. for the experimental event selection employed in 
the analyses, all the events correspond to fragments 
formed at the same temperature, and 

3. distortions on the isotopic temperature due to sec- 
ondary decay of hot primary fragments may be ne- 
glected. 

Although the Statistical Multifragmentation Model 
used in the discussion below, is based on the first as- 
sumption, the last two hypotheses are not supported by 
the model, as we shall discuss in detail. 

The SMM uses the Monte Carlo method and averages 
observables with the statistical weight over decay parti- 
tions. A multifragment decay partition is defined in the 
SMM approach || as a specific set of emitted fragments 
and light particles. For simplicity, each partition in the 
SMM approach is weighted according to the entropy of 
the partition. This entropy is approximated by analytical 
expressions rather than by an event by event sampling of 
the phase space as in ref. fl0|| These approximations rely 
upon that fact that the dominant contribution to this en- 
tropy comes from internal phase space of fragments which 
plays the role of a heat bath within the SMM approach 
just as an excited residue plays the role of a heat bath 
within compound nuclear decay theory Jl6| . 

For a given decay partition and by making a Wigner 
Seitz approximation to the Coulomb energy, energy con- 
servation within the SMM approach leads to the expres- 
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where Eq is the total excitation energy and E^ s is ground 
state energy of a nuclei having a mass and atomic num- 
ber equal to that of the total system, Aq and Zq, respec- 
tively. The first term on the right hand side stands for 
the Coulomb energy of a homogeneous charge Z^e oc- 
cupying the volume of the system of radius Rq and Naz 
indicates the number of fragments of mass number A and 
atomic number Z in the partition of the system. 

In the equation above, Eaz is the kinetic plus internal 
energy for each of these fragments. It is related to the 
temperature by assuming all fragments are at a common 
temperature as follows, 



E A 2 



±T + E* AZ (T) + E C AZ - B AZ (3) 



where the internal excitation energy of the fragments, 
E AZ (T), may be approximated by an extension of the 
semi-empirical mass formula to finite temperatures [[)), 
and the extra Coulomb energy of the fragment in the frag- 
mentation volume, E AZ , may be calculated within the 
Wigner-Seitz approximation. Baz stands for the ground 
state binding energy for the fragment. Eqs. || and || re- 
sult from an average of the microcanonical expression for 



energy conservation over the phase space corresponding 
to the specific decay partition. 

By applying the energy conservation relationship in 
Eqs. one obtains a temperature T that describes 

the internal excitation and translational energies of frag- 
ments within a given partition. Even though the overall 
system is assumed to be in equilibrium at a fixed excita- 
tion energy Eq, different decay partitions have different 
Coulomb, binding, and translational energies and, conse- 
quently, different excitation energies of the emitted frag- 
ments. Consistency with Eqs. (|^-||) therefore requires 
that the temperature T of the fragments varies from one 
decay partition to another, reflecting the differences be- 
tween the Coulomb, binding and translational energies of 
the various partitions. 

Labeling the partition {Naz} with the index, /, the 
statistical weight associated with the partition, 
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may be found by expressing the entropy of the fragments, 
Saz, using approximations derived from the liquid drop 
model at finite temperature |{^]. Consequently the phys- 
ical observables can be expressed by a weighted average 
over decay partitions as, 



(O 



AZ 



EfW f OAz 



(5) 



where O az can be any interesting observables such as the 
yield of a fragment or the temperature (In the present 
work, the summation included 10 8 partitions.) 

This allows one to predict the various results from the 
SMM that are addressed in the next section with re- 
gard to the temperature variations. Because the SMM ap- 
proach invokes a temperature to sample the microcanon- 
ical phase space, we denote the predicted observables as 
approximately microcanonical. Despite this caveat, we 
note that this procedure can, in principal, give accurate 
microcanonical predictions for experimental observables 
provided the thermal expressions for the free energies are 
accurate descriptions of the integration over the micro- 
canonical phase space. 

Before passing on to the various results of our investi- 
gation, it is important to clarify that we do not invoke the 
Grand Canonical approximations to the SMM approach 
introduced in ref. |Tjj to allow Monte Carlo event simu- 
lations Pq |. Instead, we have adhered closely to original 
SMM approach outlined in ref. (|], with the exception 
that all calculations in this paper were performed at a 
constant freezeout density equal to one third that of the 
saturation density of nuclear matter. 



2 



III. PRIMARY TEMPERATURES 

The SMM procedure expressed in Eqs.(||||) leads to a 
distribution of the temperatures of the fragmenting sys- 
tem for a given excitation energy in the same sense that 
the temperature of the daughter nucleus in compound 
nuclear decay theory varies as a function of the Coulomb 
barrier and separation energy of each decay channel. 
The points in Fig. (|l]) denote the temperature distri- 
butions for the fragmentation of an excited 112 Sn nu- 
cleus at three different excitation energies obtained with 
the SMM. These distributions are well fitted by gaus- 
sian functions, shown by the curves in the figure, with 
variances a 2 -, that are fairly independent of the energy, 
ot ~ 0.4 MeV, in the range 3MeV < E* /A < lOMeV. 
At each excitation energy, we average over all of the par- 
titions and define this average value as the " approximate 
microcanonical" temperature Tmic- 

Since each of the isotopes employed in the thermome- 
ter has a specific mass, charge and binding energy, the 
application of conservation laws sets a constraint on the 
values available to the remainder of the system. Because 
of this finite size effect, the temperature distribution ob- 
tained when a specific isotope is present is slightly differ- 
ent from the one obtained when all partitions are consid- 
ered. In particular, a small difference (< 0.1 MeV) is ob- 
served between the average temperatures for the various 
isotopes; this is illustrated in Fig. (||) for carbon isotopes 
from the fragmentation of a 112 Sn nucleus at Eq/A = 6 
MeV. Even though the average temperatures are differ- 
ent reflecting the different binding energies of the three 
isotopes, all these distributions are gaussians with nearly 
the same variances. We can extract another temperature 
Timf by averaging over partitions which contain an In- 
termediate Mass Fragment (IMF) with 3 < Z < 10. It's 
interesting to note that Tmic can exceed Timf at low 
energies by as much as 0.2 MeV, in part because it takes 
more energy to emit an IMF than to emit an equivalent 
mass in the form of alpha particles, leaving less energy 
for thermal excitation. 

The basic idea contained in Eq. (|l|) was derived under 
the assumption that the primary yields are well repre- 
sented by the grand canonical approximation at a sin- 
gle breakup temperature; the double ratio was invoked 
to cancel out the contribution to the yields coming from 
the neutron and proton chemical potentials. In the SMM, 
however, the temperature varies from partition to parti- 
tion and the chemical potentials, which appear within 
the grand canonical formalism as Lagrange multipliers 
that conserve charge and mass, are not explicitly invoked. 
Thus, we can not presume the validity of the Albergo's 
formula ( Eq. [j]) in the SMM and must test its validity 
instead. 

We begin with a test of the validity of Eq. ([l]) when 
one employs the primary yields. For a given decay parti- 
tion {Naz}, we take into account the internal free energy 
F AZ (T) which is parameterized as: 



F^z = -B(A, Z) + F* AZ {T) + F° z , 
F*az(T) = Ff|(T) + F S AZ {T) - T\n{g A s z 



(6) 
(7) 



where g AZ is the ground state spin degeneracy, and F AZ , 
F AZ , and F AZ correspond to the excitation energy de- 
pendent bulk, surface, and Coulomb contributions to the 
internal free energy |20| after the binding energy part has 
been removed. The reader is referred to ref. |9[| for ex- 
plicit expressions for the terms entering in the equation 
above. Then the primary yield for the ground state can 
be related to the total yield by 



N 9 A S Z = N AZ ■ g\ s z exp [F* AZ {T)/T] 



(8) 



for this partition. Following the procedure described in 
the previous section, we will use this expression and Eq. 
(||) to obtain the average g.s. yield distribution (N AZ ). 
This, in turn, can be used in Eq. ([!]) to extract isotopic 
temperatures as follows, 
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(9) 



In previous SMM calculations, experimental binding en- 
ergies and spin degeneracy factors g AZ were used for light 
nuclei with A < 5. Liquid drop binding energies and spin 
degeneracy factors of unity were used for A > 5. In this 
work, we will retain these conventions on spin degeneracy 
factors so as to be consistent with prior calculations, but 
we will use empirical binding energies for all nuclei. 

In Fig. (||), the isotopic temperatures Tf™ m for the 
carbon thermometer [Z\ = Z^ = 6,Ai = 11, Aa = 12) 
are plotted as the stars for the multifragmentation of a 
112 Sn source at excitation energies Eq /A = 3 — 10 MeV. 
For comparisons, the corresponding Tmic and Tjmf for 
the same system are also shown in Fig. (||) as the dashed 
and solid lines, respectively. While supporting the con- 
cept of isotopic thermometry, the good agreement be- 
tween Tjmf and Tf™ m is somewhat surprising, given the 
strong dependence of the Boltzmann factor on temper- 
ature for large AB and the width of the temperature 
distribution shown in Fig. ([!]). As shown in the follow- 
ing section, it occurs in part due to a large cancellation 
involving the Boltzmann factor and the temperature de- 
pendences of the effective chemical potentials. Fig. (||) 
also reveals that fairly precise information about Timf 
and somewhat less precise information about Tmic is 
provided by the primary yields. This suggests that given 
a precise relationship between primary to the final yields, 
it would be possible to determine the breakup tempera- 
ture from the measured yields. 



IV. EFFECTS OF TEMPERATURE VARIATIONS 

The surprising consistency between Timf and Tf™ m in 
Fig. (|3|) suggests that the corrections to the grand canon- 
ical prediction for the isotope temperatures are small, 
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and one may utilize this approach to understand why 
the temperature variations have so little influence on the 
results. Taking this tact, we assume that the isotopic 
distributions are well approximated for each partition by 
the grand canonical limit, use this limit to gain insight 
into the finite size effects and at the same time, inves- 
tigate the accuracy of this approximation. We take this 
approach to consider first the influence of the tempera- 
ture variations and later the consequences of the finite 
size on the effective chemical potentials. 

Considering the influence of the temperature varia- 
tions in this approximation, we average the grand canon- 
ical approximation over the temperature distribution in 
Fig. (|l). If the approximation works, the expressions 
that result from this average should be appropriate for 
the consideration of the effects of temperature distribu- 
tions arising from other effects and within other equilib- 
rium models of multifragmentation as well. Taking this 
approach, the yield of a particular isotope i in the frame- 
work of Albergo's method ]l5| ], when averaged over all 
possible partitions, becomes: 



V I dTf{T) 
la 



A Z J 2 UT) 



exp [(Zi fip F (T) + Ni unf (T) + B^/T] (10) 



where f(T) is the temperature distribution, V represents 
the free volume of the system, At = y // 2irh 2 /mT, m is 
the nucleon mass and /ipf (/j>nf) stands for the chem- 
ical potential associated with free protons (neutrons) at 
temperature T. The internal partition function of the 
fragment i is given by: 
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(11) 



where AEj is the excitation energy of the state j with 
respect to the ground state and g\ stands for the spin 
degeneracy factor of this excited state. 

Assuming that f(T) is a gaussian centered at (T) and 
with width op <C (T) (see Fig. [j] ), one may expand 1/T, 
T 3 / 2 , and the chemical potentials. By considering only 
fragments observed in the ground state, i.e. d(T) = gf, 
we obtain that 
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(12) 



where A* = jm. In the above expression, the cor- 

rections to the grand canonical relationship are provided 

r 2 1 

by the correction factor • exp |^ which depends on 



assumed width of the temperature distribution, the bind- 
ing energy of the i-th fragment, the neutron and proton 
chemical potentials and their derivatives through the pa- 
rameters p and q. These two parameters are defined by 
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The isotopic temperature can be extracted from the 
above corrected yields. Replacing Y(A, Z) in Eq. (|l|) by 
the right hand side of Eq. (|l2|), one cancels out the spin 
and mass dependent term C and then obtains: 



exp [AB/T%J] 
where 

G(A, Z) = exp 
1 



G{A U Z 1 )/G{A 1 + \,Z 1 ) 
G(A 2 ,Z 2 )/G(A 2 + 1,Z 2 ) 



(15) 
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• exp 
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(16) 



In the above double ratio the terms involving the chemi- 
cal potentials evaluated at the average temperature can- 
cel; however, terms in the correction factor involving the 
derivatives of the chemical potentials remain. 

Quantitative estimates of the correction factor require 
one to obtain estimates for the effective chemical poten- 
tials and their derivatives with respect to temperature. 
The proton and neutron chemical potentials at tempera- 
ture T may be calculated from the free proton and neu- 
tron multiplicities via the expression: 



Hpf(T) = TTog 



(J-nf{T) = Tlog 



r A T Y PF (T) 



9pfV 
r\ 3 T Y NF (T) 



(17) 



9nfV 



where gpF(gNF) represents the spin degeneracy factor 
of the proton(neutron). For the calculations presented 
in this work, it has proven advantageous and reasonably 
accurate to approximate the yields Ypp{T) and Ynf(T) 
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over a modest range in temperature by power law expres- 
sions in the temperature. In this approximation, 



Y PF (T) = cpfTI" 
Y NF (T) = c NF T^ 



(18) 



For the decay of 112 Sn nuclei at temperatures ranging 
over 4 < T < 7 MeV , Yp F and Y/vf are well described 
by 7pf = 4.5 and jnf = 1-0 {cpf = 1-33 x 10~ 4 and 
cnf = 0.267) according to the SMM; comparisons of 
this parameterization to yields calculated with the SMM 
model are shown Fig. (Q). These values depend on the 
density, which has been chosen to be one third that of 
the saturation density of nuclear matter. Larger values 
of the free nucleon yields are obtained at lower density. 

Using this approximation, the explicit forms of the cor- 
rection factors in Eqs. (|l^-([l4|) become 2a FF — f3p F = 
(iPF - |) = 3 and 2a NF = f3 NF = (~f NF - §) = -|. 
We note that the correction factor to the temperature 
in Eq. (|l|) depends on the power law exponents 
1pf{inf) and not on the absolute values of the pro- 
ton(neutron) yields. 

Even though Eq. ( |l0| ) has an exponent that appears 
to be strongly temperature dependent, there is a strong 
cancelation between the contributions from the chemical 
potentials and binding energy factors in the expressions 
for p and q. As a result, the correction factor is of order 



unity. Values in the range of -7= • exp 



1 — 2 are 



obtained, for example, in the decay of Sn nuclei at 
temperatures in the range of 4 < T < 7MeV . 

The isotopic temperatures T?£f calculated from Eq. 
( IHf ) for carbon thermometer are shown in Fig. (||) in 
comparisons with temperatures Tmic, Timf and T?™ m 
derived from the SMM in the previous session. The very 
good agreement between T^J , T?™ m and T IM f indi- 
cates that the corrections to the isotopic temperatures 
associated with these temperature variations are small, 
although the yields can change by as much as a factor 
of two. This comparative insensitivity arises because 
the isotopic thermometers depend logarithmically on the 
yields. 

This insensitivity depends on the nature and magni- 
tude of the temperature variation. The corrections to the 
isotopic temperatures will be somewhat larger in other 
contexts or other models where the temperature varia- 
tions are larger. The limited precision with which sys- 
tems may be selected experimentally may also have a 
similar influence because the excitation energy and tem- 
perature varies experimentally from collision to collision 
due to variations in the impact parameter or in the en- 
ergy removed by preequilibrium particle emission. The 
influence of this temperature variation, which may ex- 
ceed the variation in temperature caused by the aver- 
aging over decay partitions, can also be estimated via 
techniques outlined in the present section. To illustrate 
how one can estimate the possible corrections due to an 
imprecision in the excitation energy definition, the circles 



in Fig. (g) show calculations using Eq. ( jT5| ) for carbon 
thermometer assuming a width of a F ~ 0.8 MeV for the 
temperature distribution, which is twice as large as that 
predicted in Figs. 1 and 2. This width is not based upon 
a dynamical calculation; it is only to illustrate that larger 
isotopic temperatures can result if the excitation energy 
is poorly defined. 



V. CHEMICAL POTENTIALS 

The grand canonical limit has a great advantage of 
providing a simple analytic expression for the isotopic 
yields from which other useful expressions can be de- 
rived. However, the concept of uniform chemical poten- 
tials is not a prediction of microcanonical or canonical 
models and must be investigated to determine its appli- 
cability to finite systems. We do this by trying to com- 
paring the grand canonical expression for the isotopic 
yields to the predictions of approximately microcanoni- 
cal SMM calculations. We start by assuming that these 
isotopic distributions can be calculated within the grand 
canonical approximation and then test this assumption 
as follows. Using a pair of adjacent isotopes, we invert the 
grand canonical expression for the isotopic yields of two 
adjacent isotopes to obtain an equation for the effective 
neutron chemical potential: 



tf f (A,Z) =Tlog[ 



as 

9az 



9 g A + iz U+l 



.4 



3/2 



exp((B AZ -B A+lz )/T) 



1 AZ 



(19) 



where g AZ , B A z and YjfL are the ground state spin de- 
generacy, the binding energy and the ground state pri- 
mary yield for a fragment with (A,Z), respectively. If the 
Y^ z taken to be the ground state yields predicted by the 
SMM, 11° ff (A, Z) becomes an effective "SMM" chemical 
potential. By performing SMM calculations, we find the 
temperature- and isotopic- dependences of the effective 
neutron chemical potentials given in Fig. (|) for Carbon 
and Lithium isotopes from the decay of a Sn nucleus 
at excitation energies of Eq j A = 3, 6, 9 MeV. 

These effective chemical potentials are essentially the 
same for the Carbon and Lithium isotope chains. This 
insensitivity to element number is consistent with the 
concept of a chemical potential and offers support for 
the use of the grand canonical expression to describe iso- 
topic distributions. There is a dependence on the neu- 
tron number of the isotope, however, that lies outside 
of the grand canonical approximation. This variation in 
the neutron chemical potential basically comes as a result 
of mass, charge and energy conservation for a finite-size 
system. We can understand the influence of these conser- 
vation laws most easily at low excitation energies, where 
the two largest fragments in the final state are the IMF 
(Carbon or Lithium in this case) and a heavy residue 
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which contains most of the remaining charge and mass. 
We estimate the influence of conservation laws at low ex- 
citation energy qualitatively by considering binary decay 
configurations. Assuming that a parent nucleus (Aq, Zq) 
decays into a light fragment (A,Z) and a heavy residue 
(Aq — A, Zq — Z) , we can approximate the yield of frag- 
ment (A,Z) in its ground state by 



Y&ctpP' (A,Z)p* (A 



A,Z -Z)p REL (20) 



g g A s z -exp[S* (A -A,Z -Z)} 



A-(A - A) 



A 



3/2 



A-(Aq-A) 



where p gs — g AZ , p* and S* are the density of states 
for the light nucleus in its ground state level, the den- 
sity of states and entropy of the heavy residue in its 
excited state, respectively, The other factor, /Orex, ~ 

3/2 

X T , is the thermal average of the state 

density of relative motion. 

Replacing the yields in Eq.([l9|) with Eq.(p0|) and as- 
suming A << A , one finds that the effective chemical 
potentiail depends on the difference in residue entropies, 
S* (Ao - A - 1, Z) - S* (A -A,Z -Z). Using an ex- 
pansion for small changes in the nuclear entropy from 
ref. (T^| , this difference can be expressed in terms of the 
difference of binding energies, 

S* (A - A- 1, Z) - S* (A -A,Z - Z) 
= -{Ba -a,z -z - Ba -a-\,z -z)/T 

-{B A z-B A+ iz)/T + f*/T (21) 

plus a term depending on the free excitation energy per 
nucleon, /* = E* /A - TS/A . This difference in bind- 
ing energies is further related to the neutron separation 
energy s n (A - A,Z - Z): 

s n {Ao — A, Z — Z) = BA a -A.Za-z — Ba -a-i,Zo-z 

(22) 

One consequently obtains the following expression for the 
effective chemical potential: 



Pn 



-s n (A -A,Zo-Z)+f* 



(23) 



where the reduced free excitation energy has been ap- 
proximated by its low energy limit, 



/* = - — , £ = 8MeV. 



(24) 



12 C+X , the chemical potential 
A, Zq — Z), is plotted as the 



For the decay 112 Sn — i 
at T = 0, i.e., — s n (Ao 

stars in Pig. (pi) ; the binding energies for these calcula- 
tions were calculated using the liquid-drop parametriza- 
tion in ref. [ pl| . The reduced free energy /* gives a rea- 
sonable estimate for the trend with excitation energy 



The dot-dashed line in Fig. (Bh gives the chemical po- 
tential predicted from Eq. (||) for Eg /A = 3 MeV 
(T = 4.58MeV). The predicted trend is close to that 
predicted by the SMM model (solid circles and squares) 
but has a somewhat stronger dependence on N — Z. 

In general, the slope of the effective neutron chemical 
potential is getting slightly flatter as the excitation en- 
ergy or temperature increases. If we consider that the 
system undergoes a multiple fragment decay at higher 
temperatures, it is clear that approximating the entropy 
of the remaining system by that of a residue of compa- 
rable mass becomes rather inaccurate. The constraints 
imposed on the total system by the isospin asymmetry of 
one observed fragment should, in that case, be less signif- 
icant. While there is a mass dependence to the effective 
chemical potential that is inconsistent with the grand 
canonical approach, it is useful to note that the mass de- 
pendence of the chemical potential (for these systems of 
more than 100 nucleons) is small if one is mainly con- 
cerned with nuclei near the valley of stability. If one can- 
cels the chemical potential effects by constructing double 
ratios like that of the Algergo formula, the consequence 
of such finite size effects becomes negligible indeed. 



VI. INFLUENCE OF SECONDARY DECAY 



As discussed in sect. D, fragments are formed in ex- 
cited states as well as in their ground states, correspond- 
ing to the breakup temperature. Fragments in short lived 
excited states decay before they are detected and, there- 
fore, the observed yields differ from that of the primary 
fragments. The effects of secondary decay on the iso- 
topic yields and isotopic temperatures have already been 
reported by some authors (see for example [p2|-p4|). Al- 
though the approaches employed in the description of the 
decay of hot primary fragments are different, all those 
works qualitatively agree on the point that the isotopic 
temperature is lower than the thermodynamical one. 

At the quantitative level, details of the population and 
decay of the excited fragments are important. One issue 
concerns the importance of utilizing empirical binding 
energies, energy levels and decay branching ratios for the 
excited fragments. Fig. (^J) shows the primary and sec- 
ondary carbon isotopic distributions for the decay of a 
112 Sn nucleus at initial excitation energies of Eq/A = 4 
and 6 MeV. The primary distribution (solid line) is cal- 
culated by considering empirical binding energies within 
the SMM for hot fragments. The si mp lified Weisskopf 
evaporative decay procedure of ref. []17| is used for one 
final distribution (dotted line). The other final distri- 
bution (dashed line) is obtained by calculating the sec- 
ondary decay for Z < 10 hot fragments, as in ref. p4] , p5[ , 
according to empirical nuclear structure information re- 
garding the excitation energies, spins, isospins and decay 
branching ratios where available. For hot fragments with 



G 



Z < 10 where such information is not available, the de- 
cay is calculated according to the Hauser-Feshbach for- 
malism |26|]. The contributions to this latter calculation 
from the secondary decay of hot fragments with Z > 10, 
are calculated, for simplicity, via the secondary evapo- 
rative decay procedure of ref. [0. Decays of fragments 
with Z > 10 make a 15% contribution to the yields of 
12 C isotopes that may be altered when the decay of hot 
fragments with Z > 10 is calculated more accurately. 

Obviously, in Fig. (||), the final distribution after the 
empirical secondary decay is much wider than the final 
distribution obtained via the evaporative decay approach 
of ref. |17j]. This points out the importance of using 
the empirical information in such calculations. This also 
leads to the extraction of larger isotopic temperatures via 
Eq. (^) for the empirical approach. Temperatures for 
the Carbon isotope thermometer and He-Li thermome- 
ter calculated for the two secondary decay approaches 
are shown, for example, in Fig. ([7|) for the multifrag- 
mentation of a 112 Sn nucleus at Efi/A = 4 - 10 MeV. 
For reference, the curves Tmic an d Timf from Fig. (||) 
are also shown as the dashed and solid lines in the figure. 
Clearly, incorporating empirical information in the decay 
makes a significant difference. Both calculations provide 
lower isotopic temperatures than have been obtained in 
recent experiments [fj|||, |l3| , [L4] . 

It should be noted, however, that the simplified Weis- 
skopf evaporative decay, shown in Figs. (||) and (|7|), is 
only used in the SMM code of ref. |l7]] to calculate the de- 
cay of fragments with A > 16. The decay of lighter frag- 
ments is calculated via a "Fermi Breakup" multiparticle 
decay formalism. This latter decay mechanism makes 
the dominant contribution to the isotope temperatures 
calculated via the latter SMM code in ref. |27j . Investi- 
gations of the experimental and theoretical basis for the 
" Fermi Breakup" approach are needed, but are out of the 
scope of the present work. 

Regardless of the decay formalism, memory of the 
breakup stage is lost via the secondary decay mechanism. 
The degree of memory loss depends on the details of the 
secondary decay correction and on the role of short-lived 
higher lying particle unbound states. A smaller degree 
of memory loss ensues in models such as those of refs. 
jl0|,^8 129|, where few, if any, particle unbound states are 
considered. The approach of ref. 1 17| represents the other 
extreme, wherein all states are considered regardless of 
lifetime. This issue clearly needs further study to see 
whether the role of particle unstable nuclei can be con- 
strained, for example, by direct measurements using tech- 
niques discussed in refs. J|5|,|30| or by other experimental 
observables. 



VII. CONCLUDING REMARKS 



We discussed some of main aspects that could cause 
microcanonical predictions for isotopic distributions and 



isotopic temperatures to differ from grand canonical cal- 
culations and influence the determination of the breakup 
temperature and other experimental observables. We in- 
vestigate this problem by checking the consistency of the 
grand canonical expression for the isotopic yields against 
the approximately microcanonical SMM predictions and 
explore the potential role which may be played by vari- 
ations in the temperature and in the effective chemical 
potentials. These variations occur as a consequence of 
the finite size of the disintegrating system and are there- 
fore present in all microcanonical calculations. 

Concerning the temperature variation, we find that 
this causes the isotopic yields obtained with the approxi- 
mately microcanonical SMM simulations for the primary 
distribution to differ from those of the grand-canonical 
ensemble by factors of order unity. One difference stems 
from the averaging over the temperatures corresponding 
to the different breakup partitions. These vary because 
the total binding, coulomb and translational kinetic ener- 
gies vary from partition to partition and by subtraction, 
the thermal energy must vary as well. A simple and rela- 
tively accurate prescription that accounts for these tem- 
perature variations was given that may also prove useful 
for estimating the influence of thermal averaging over 
the variations in the actual excitation energy deposition 
within a data set that is constrained by an experimental 
cut on the estimated energy deposition. 

We also extract effective chemical potentials by com- 
paring approximate microcanonical and grand canonical 
expressions for the isotopic yields. These effective chem- 
ical potentials are approximately the same for isotopes 
of different elements that lie along the valley of beta sta- 
bility, but vary as a function of (N-Z). For example, we 
observe for the neutron chemical potential a dependence 
upon (N-Z) that can be understood at low excitation en- 
ergies to arise from the dependence of the neutron separa- 
tion energy on the location of the accompanying residue 
relative to the line of beta stability. 

Typically, these variations in temperature and effective 
chemical potential cause variations in the isotopic yields 
of order unity. The logarithmic relation between the iso- 
topic temperature and the yields means that the latter 
may be wrongly predicted by a factor of two and one 
may still find a reasonable agreement between the ap- 
proximate microcanonical and the isotopic temperatures 
provided the binding energy difference AB is significantly 
larger than the temperature. When the effects of sec- 
ondary decay is taken into account, however, the yields 
can change by more than an order of magnitude and the 
temperature values can decrease appreciably. While the 
magnitude of this change is not yet unambiguously estab- 
lished, it was shown that the incorporation of empirical 
information about the decay is essential for quantitative 
comparisons to experimental data. Measurements that 
quantify the role of higher lying particle unstable states 
are essential for determining the magnitude of these sec- 
ondary decay corrections. 
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FIG. 1. The points denote distributions of temperatures 
calculated with the SMM approach for the decay of a 112 Sn 
nucleus at three different excitation energies. The lines denote 
gaussian fits to the calculated distributions. 

FIG. 2. The points denote temperature distributions calcu- 
lated with the SMM approach for the different isotopes con- 
sidered in the carbon thermometer for an excitation energy 
of Eq I A — 6MeV. The lines denote gaussian fits to the calcu- 
lated distributions. 



FIG. 3. Comparisons of various primary temperatures 
T MIC , Timf and T?£ m from the SMM and from the an- 
alytical calculation in the grand canonical limit. For details 
see the text. One point is missing for with ot = 0.8MeV 
because the calculated value for p for the correction term in 
Eq.(p^) becomes negative at Eq/A — 3MeV, i.e. the expan- 
sion breaks down in this case. 



FIG. 4. The solid squares and circles denote the free pro- 
ton and neutron yields, respectively, calculated via the SMM 
approach. The solid and dashed lines denote fits to the cal- 
culated yields following Eq. (|l8|). 

FIG. 5. The squares, circles and triangles denote neutron 
chemical potentials derived from Eq. (y_9|) using SMM pre- 
dictions for Carbon and Lithium isotopic yields at various 
initial excitation energies for the decay of the nucleus 112 Sn. 
The stars and the dot-dashed line denote approximate values 
calculated from Eq.Q for T=0 and 4.58 MeV, respectively. 
The error bars denote the statistical errors in the calculation, 
which in many cases are too small to be observed in the figure. 

FIG. 6. Primary (solid line) and final Carbon isotopic dis- 
tributions calculated for the decay of the nucleus 112 Sn using 
(dashed line) and neglecting (dotted line) the empirical nu- 
clear structure information in the secondary decay process. 
The error bars denote the statistical errors in the calculation, 
which in many cases are too small to be observed in the figure. 

FIG. 7. Isotopic temperatures for Carbon and He-Li ther- 
mometers calculated with the SMM model for the decay of 
the nucleus 112 Sn using (solid symbols) and neglecting (open 
symbols) the empirical nuclear structure information in the 
secondary decay process. The lines are the same as those 
shown in Fig. (ph. The error bars denote the statistical errors 
in the calculation, which in many cases are too small to be 
observed in the figure. 
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